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Abstract 

An algorithm for an improved description of final-state QCD radiation is 
introduced. It is matched to the first-order matrix elements for gluon emis- 
sion in a host of decays, for processes within the Standard Model and the 
Minimal Supersymmetric extension thereof. 



The objective of this article is to summarize the improved description of parton shower 
evolution Jl| recently introduced starting with Pythia 6.154 ||. In particular, process- 
specific 0(a s ) matrix elements for gluon emission in decays a —> be are used to match the 
shower description to the correct emission rate in the hard-gluon region, and to provide 
the proper amount of 'dead cone' || suppression of collinear gluon emission off massive 
particles. The original motivation was to improve the understanding of bb events at LEP1. 
For linear colliders the applications to top, Higgs and SUSY physics are very important. 

The traditional final-state shower algorithm [|J in Pythia is based on an evolution 
in Q 2 = m 2 , i.e. potential branchings are considered in order of decreasing mass. A 
branching d — > ef is then characterized by m\ and z = E e /Ed- For the process 7*/^ — ► qq, 
the first gluon emission off both q and q are corrected to the first-order matrix elements 
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for 'y* /Z — > ggg. (The a s and the Sudakov form factor are omitted from the comparison, 
since the shower procedure here attempts to include higher-order effects absent in the 
first-order matrix elements.) 

This matching is well-defined for massless quarks, and was originally used unchanged 
for massive ones. A first attempt to include massive matrix elements did not compensate 
for mass effects in the shower kinematics, and therefore came to exaggerate the suppression 
of radiation off heavy quarks || . Now the shower has been modified to solve this issue, 
and also improved and extended better to cover a host of different reactions [p]]. 

The starting point is the calculation of processes a — > be and a — > beg, where the ratio 

txt i \ 1 da(a ^ beg) 

W ME (x l ,x 2 ) = - — \ 1 

a{a — > be) dx\ dx 2 

gives the process-dependent differential gluon-emission rate. Here the phase space vari- 
ables are x% = 2Eb/m a and x 2 = 2E c /m a , expressed in the rest frame of parton a. Using 
the standard model and the minimal supersymmetric extension thereof as templates, a 
wide selection of colour and spin structures have been addressed, exemplified by Z° —> qq, 
t -> bW + , H° -» qq,t -»■ bH + , Z° -> qq, q -> q'W + , H° -> qq, q -> q'H+, % -> gg, 
g — > qx, t ^ ix, g —> qq, q ^> qg, and t — > tg. The mass ratios 77 = m b /m a and 
r 2 = m c /m a have been kept as free parameters. When allowed, processes have been cal- 
culated for an arbitrary mixture of "parities", i.e. without or with a 75 factor, like in the 
vector/axial vector structure of J*/Z. All the matrix elements are encoded in the new 
function PYMAEL (NI , XI , X2 , R1 , R2 , ALPHA) , where NI distinguishes the matrix elements 
and ALPHA is related to the 75 admixture. 

In order to match to the singularity structure of the massive matrix elements, the 
evolution variable Q 2 is changed from m 2 to m 2 — rn 2 n _ shell , i.e. 1/Q 2 is the propagator 
of a massive particle. Furthermore, the z variable of a branching needs to be redefined, 
which is achieved by reducing the three-momenta of the daughters in the rest frame of 
the mother. For the shower history b — > bg this gives a differential probability 

«/ / , a s dQ 2 2dz 1 a s 2 

W PSA ( Xl , x 2 ) = ^C F — — — - = -C F X3{1 + r2 _ r2 _ X2) , (2) 

where the numerator 1 + z 2 of the splitting kernel for q — > qg has been replaced by a 
2 in the shower algorithm. For a process with only one radiating parton in the final 
state, such as t — > bW + , the ratio Wme/Wps,i gives the acceptance probability for an 
emission in the shower. The singularity structure exactly agrees between ME and PS, 
giving a well-behaved ratio always below unity. If both b and c can radiate, there is a 
second possible shower history that has to be considered. The matrix element is here 
split in two parts, one arbitrarily associated with b — > bg branchings and the other with 
c — > eg ones. A convenient choice is Wme.i = Wme(1 + r\ — r| — x\)/xz and Wme.2 = 
W / me(1+^1 — r i — x 2)/x^, which again gives matching singularity structures in Wme^/Wps,* 
and thus a well-behaved Monte Carlo procedure. 

Also subsequent emissions of gluons off the primary particles are corrected to Wme- 
To this end, a reduced-energy system is constructed, which retains the kinematics of the 
branching under consideration but omits the gluons already emitted, so that an effective 
three-body shower state can be mapped to an {x\, x 2 , r x , r 2 ) set of variables. For light 
quarks this procedure is almost equivalent with the original one of using the simple uni- 
versal splitting kernels after the first branching. For heavy quarks it offers an improved 
modelling of mass effects also in the collinear region. 
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Figure 1: Gluon radiation pattern Wme(%i, x 2), eq. (|l|), as a function of the gluon emission 
angle, for four fixed gluon energy fractions, x 3 = 2E g /E CM . The daughters/mother mass 
ratio is fixed at 0.2. The curves are for different combinations of colour and spin, and 
with/without a 75 factor where allowed. 



Some further changes have been introduced, a few minor as default and some more 
significant ones as non-default options [|TJ. This includes the description of coherence 
effects and a s arguments, in general and more specifically for secondary heavy flavour 
production by gluon splittings. 

Further issues remain to be addressed, e.g. radiation off particles with non-negligible 
width. In general, however, the new shower should allow an improved description of gluon 
radiation in many different processes. Where it can be tested, for the amount of radiation 
off b quarks relative to light ones at LEP1, the new algorithm indeed is successful [[[J [|. 

As an illustration of the process dependence, Fig. [I] shows the radiation pattern of 
the various matrix elements calculated. In order to ease the comparison, the same fixed 
mass ratios have been used for all processes, r\ = r 2 = 0.2. Furthermore, the large mass 
ratio highlights the dead cone effect, which shows a universal behaviour for small gluon 
energies. At large angles, and still small gluon energies, there is a dependence on the 
colour structure of the process, but not e.g. on the spin of the particles. This should be 
expected, since in the soft-gluon limit radiation can be described by a spin-independent 
eikonal expression |J. Maybe more surprising is how completely this universality breaks 
down for more energetic gluons. Then processes are split not only by colour, but also by 
the spin structure, and the presence or not of a 75 in the matrix element, where allowed. 
(The figure only show the two extremes; by an arbitrary admixture of the two one would 
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instead obtain a set of allowed bands.) Furthermore, the dead cone effect is shown to 
remain only for the case of a spin particle decaying to two daughters also with spin 0. In 
retrospect, the process dependence is there also at small gluon energies, but is nonsingular 
and therefore invisible underneath the eikonal soft-gluon-singular contributions. 

The above figure well illustrates that differences could be big in principle, but fortu- 
nately the reality is more forgiving. One reason is the big jump in mass between the b 
quark, on the one hand, and t, SUSY and any other potentially coloured particles, on the 
other. The most direct consequence is that the heavier particles typically generate only a 
small fraction of the total amount of QCD radiation, while b and lighter quarks produce 
the bulk of it. The b is light on the scale of the decaying particle, and so has a smaller 
dead cone than the one in Fig. [I]. 

A more realistic example of differences is then offered by a light Higgs state, say 115- 
130 GeV in mass as suggested by the MSSM scenario, decaying to bb. The three-jet rate in 
such events typically is 10-15% higher than in 7*/Z* — > bb (or light quark) decays at the 
same energy. The difference is less for soft radiation, so the Higgs decay is only producing 
about 1% lower mean values for the b quark and B hadron fragmentation functions. 

In tt events, the new algorithm increases the amount of radiation in the top production 
stage, but decreases it in the subsequent top decay. The difference is especially notable 
in the W hemisphere of the top decay, where the gluon emission rate is dropping rather 
steeper (with the angle away from the b quark) in the new program than in the old. This 
is related to a destructive interference between emission off the t and off the b in this 
hemisphere, while the older approach had its origins in e + e~ — > qq events, where the 
interference is constructive. The net result is a small but visible decrease in the total 
amount of gluon radiation in tt events. 

For supersymmetric processes, results largely depend on the actual masses. Assuming 
the charged Higgs mass to coincide with the one, the decays t —> bW + and t — > bH + 
give almost identical amounts of radiation. But if the stop mass agrees with the top one, 
there is more QCD radiation in the former production process than the latter. (While the 
difference in threshold behaviour here gives the opposite effect for ISR photon radiation, 
which can become more important.) 
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